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Abstract 

We study the effective field theory of inflation, i. e. the most general theory describing the fluctuations 
around a quasi de Sitter background, in the case of single field models. The scalar mode can be eaten 
by the metric by going to unitary gauge. In this gauge, the most general theory is built with the 
lowest dimension operators invariant under spatial diffeomorphisms, like g^^ and Kfj_^, the extrinsic 
curvature of constant time surfaces. This approach allows us to characterize all the possible high 
energy corrections to simple slow-roll inflation, whose sizes are constrained by experiments. Also, it 
describes in a common language all single field models, including those with a small speed of sound 
and Ghost Inflation, and it makes explicit the implications of having a quasi de Sitter background. 
The non-linear realization of time diffeomorphisms forces correlation among different observables, 
like a reduced speed of sound and an enhanced level of non-Gaussianity. 



1 Introduction 



The effective held theory approach, i.e. the description of a system through the lowest dimension 
operators compatible with the underlying symmetries, has been very fruitful in many areas, from 
particle physics to condensed matter. The purpose of this paper is to apply this methodology to 
describe the theory of fluctuations around an inflating cosmological background. 

The usual way to study a single field infiationary model is to start from a Lagrangian for a scalar 
field 4> and solve the equation of motion for (j) together with the Priedmann equations for the FRW 
metric. We are interested in an inflating solution, i.e. an accelerated expansion with a slowly varying 
Hubble parameter, with the scalar following an homogeneous time-dependent solution (j)o{t). At this 
point one studies perturbations around this background solution to work out the predictions for the 
various cosmological observables. 
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The theory of perturbations around the time evolving solution is quite different from the theory 
of (j) we started with: while (/> is a scalar under all diffeomorphisms (diffs), the perturbation 5(j) is a 
scalar only under spatial diffs while it transforms non-linearly with respect to time diffs: 

+ 5<l) ^ 5<l) + ^Q{t)i^ . (1) 

In particular one can choose a gauge 4>{t,x) = 4>o{t) where there are no inflaton perturbations, but 
all degrees of freedom are in the metric. The scalar variable 5(p has been eaten by the graviton, which 
has now three degrees of freedom: the scalar mode and the two tensor helicities. This phenomenon 
is analogous to what happens in a spontaneously broken gauge theory. A Goldstone mode, which 
transforms non-linearly under the gauge symmetry, can be eaten by the gauge boson (unitary gauge) 
to give a massive spin 1 particle. The non-linear sigma model of the Goldstone can be embedded 
and UV completed into a linear representation of the gauge symmetry like in the Higgs sector of 
the Standard Model. This is analogous to the standard formulation of inflation, where we start 
from a Lagrangian for (j) with a linear representation of diffs. In this paper we want to stress the 
alternative point of view, describing the theory of perturbations during inflation directly around 
the time evolving vacuum where time diffs are non-linearly realized. This formalism has been firstly 
introduced, for a generic FRW background, in [1] to study the possibility of violating the Null Energy 
Condition; here we will extend this formalism focusing on an inflationary solution. 

We will show that in unitary gauge the most generic Lagrangian with broken time diffeomor- 
phisms (but unbroken spatial diffs) describing perturbations around a flat FRW with Hubble rate 
H{t) is given by 

S = ji>x^\^-Ml^R + Ml,Hg«'-Ml^{3H^ + H)+^i^(r + lf (2) 

The first two operators after the Einstein-Hilbert term are fixed by the requirement of having a given 
unperturbed solution H(t), while all the others are free and parametrize all the possible different 
theories of perturbations with the same background solution. As time diffs are broken one is allowed 
to write any term that respects spatial diffs, including for example g^^ and the extrinsic curvature 
K^i, of the surfaces at constant time. The coefficients of the operators will be in general time 
dependent. The reader may be worried by the use of a Lagrangian that is not invariant under 
diffeomorphisms. But clearly diff. invariance can be restored as in a standard gauge theory. One 
performs a time-diffeomorphism with parameter (,^{t, x) and promotes the parameter to a field 7r(t, x) 
which shifts under time diffs: 7r(t, x) — > 7r(t, x) — C^{t, x). The scalar tt is the Goldstone mode which 
non linearly realizes the time diffs and it describes the scalar perturbations around the FRW solution. 

It is well known that the physics of the longitudinal components of massive gauge bosons can be 
studied, at sufficiently high energy, concentrating on the scalar Goldstone mode (equivalence theo- 
rem). The same is true in our case: for sufficiently high energy the mixing with gravity is irrelevant 
and we can concentrate on the Goldstone mode. In this regime the physics is very transparent and 
most of the information about cosmological perturbations can be obtained. Performing the broken 
diff transformation on the Lagrangian ([2|) and concentrating on the Goldstone mode vr one gets 



Stt= d X 



(3) 
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Every invariant operator in unitary gauge is promoted to a (non-linear) operator for the Gold- 
stone: the non-linear realization of diff invariance forces the relation among various terms. 

Let us briefly point out what are the advantages of this approach before moving to a systematic 
construction of the theory. 

• Starting from a "vanilla" scenario of inflation with a scalar fleld with minimal kinetic term 
and slow-roll potential, we have parameterized our ignorance about all the possible high en- 
ergy effects in terms of the leading invariant operators. Experiments will put bounds on the 
various operators, for example with measurements of the non-Gaussianity of perturbations 
and studying the deviation from the consistency relation for the gravitational wave tilt. In 
some sense this is similar to what one does in particle physics, where one puts constraints on 
the size of the operators that describe deviations from the Standard Model and thus encode 
the effect of new physics. 

• It is explicit what is forced by the symmetries and by the requirement of an inflating back- 
ground and what is free. For example eq. ([3]) shows that the spatial kinetic term (Vtt)'^ is 
proportional to H, while the time kinetic term tt^ is free. Another example is the unitary 
gauge operator (^f'^^ -|- 1)^. Once written in terms of the Goldstone vr, this gives a quadratic 
term vr^, which reduces the speed of sound of tt excitations, and a cubic term 7r(V7r)^, which 
increases the interaction among modes, i.e. the non-Gaussianity. Therefore, barring cancella- 
tions with other operators, a reduced speed of sound is related by symmetry to an enhanced 
non-Gaussianity. Notice moreover that the coefficient of this operator is constrained to be 
positive, to avoid propagation of vr excitations out of the lightcone. 

• One knows all the possible operators. For example, at the leading order in derivatives, the 
interaction among three vr modes can be changed by {g^^ + 1)^ and {g^^ + 1)^. This will corre- 
spond to two different shapes of the 3-point function which can be in principle experimentally 
distinguished to fix the size of each operator. 

• All the possible single field models are now unified. For example there has been interest 
in models with a modified Lagrangian L((50)^, (/>), like DBI infiation [21 [3l SJ EJ [6] which 
have rather peculiar predictions. In our language these correspond to the case in which the 
operators {g^^ + 1)" are large. Another interesting limit is when i7 ^ 0; in this case the 
leading spatial kinetic term is coming from the operator proportional to and it is of the 
form (V^vr)^. This limit describes Ghost Infiation [7]. 

• In the (j) language one can perform a field redefinition cf) — > (/>((/)). It is true that the resulting 
Lagrangian will describe the same physics, but this is not obvious. A simple example is given 
by the Lagrangian 

f{ct>)\dct>f-V{ct>), (4) 

where / is a generic function. This is equivalent to a Lagrangian with minimal kinetic term 
and a different potential through the field redefinition </)(^), d4>/d4> = f{4>). However the 
equivalence among different Lagrangians becomes more complicated when we consider more 
general terms. On the other hand this ambiguity is absent at the level of vr, which realizes a 
sort of standard non-linear representation of time diffs. 
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• In the (j) language is it not obvious how to assess the importance of an operator for the study 
of perturbations, because some of the legs of an operator may be evaluated on the background 
solution. For example in a theory with all operators of the form (c?(/))^", all of them may 
have the same importance if the background velocity <j)Q is large enough, as it happens in 
DBI inflation. On the other hand the usual way of estimating the importance of an operator 
works in the vr language. Even more clear is the case of Ghost Inflation where, given the non- 
relativistic dispersion relation for vr the scaling of operators is clear only in the vr language. 

• The parametrization of the operators directly around the solution is crucial if one calculates 
loop corrections of cosmological perturbations. A diagram with a given number of external 
legs will in general contain a UV divergence. This is easy to renormalize in the Lagrangian 
([2|), because there is only a finite number of terms which describe the interaction among n 
perturbations. On the other hand at the level of the (j) Lagrangian, there is an infinite number 
of operators contributing to the interaction among n perturbations. For each operator in fact 
one can put many of its legs on the background, so that the relation among an operator and 
a diagram for perturbations is rather obscure. 

2 Construction of the action in unitary gauge 

Inflation is a period of accelerated cosmic expansion with an approximately constant Hubble pa- 
rameter. This quasi de Sitter background has a privileged spatial slicing, given by a physical clock 
which allows to smoothly connect to a decelerated hot Big Bang evolution. The slicing is usually 
realized by a time evolving scalar (j){t). Another example one may keep in mind is given by a perfect 
fluid 0. To describe perturbations around this solution one can choose a gauge where the privileged 
slicing coincides with surfaces of constant i, i.e. 6(j){x,t) = 0. In this gauge there are no explicit 
scalar perturbations, but only metric fluctuations. As time diffeomorphisms have been fixed and 
are not a gauge symmetry anymore, the graviton now describes three degrees of freedom: the scalar 
perturbation has been eaten by the metric. 

What is the most general Lagrangian in this gauge? One must write down operators that are 
functions of the metric g^u, and that are invariant under the (linearly realized) time dependent 
spatial diffeomorphisms — > + x). Spatial diffeomorphisms are in fact unbroken. Besides 
the usual terms with the Riemann tensor, which are invariant under all diffs, many extra terms 
are now allowed, because of the reduced symmetry of the system. They describe the additional 
degree of freedom eaten by the graviton. For example it is easy to realize that g^^ is a scalar under 
spatial diffs, so that it can appear freely in the unitary gauge Lagrangian. Polynomials of g^^ are 
the only terms without derivatives. Given that there is a preferred slicing of the spacetime, one is 
also allowed to write geometric objects describing this slicing. For instance the extrinsic curvature 
K^i, of surfaces at constant time is a tensor under spatial diffs and it can be used in the action. 
Notice that generic functions of time can multiply any term in the action. In appendix |A] we prove 

"'^ Indeed, as shown for example in [8], non-vorticous excitations of a perfect fluid may be described by a 
derivatively coupled scalar. 
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that the most generic Lagrangian can be written as 

S = Jd^x ^[^Ml,R - <;(()9" - A(() + ^hhmg'" + 1)' + ^hhit fig'" + i f + 

^Mlf^g«> + _ 'JlfiK^^ - ^SK^.SK", + ...] . (5) 

where the dots stand for terms which are of higher order in the fluctuations or with more derivatives. 
We denote by SK^i, the variation of the extrinsic curvature of constant time surfaces with respect to 
the unperturbed FRW: dK^^, = K^^, — a^Hh^y with /i^i, is the induced spatial metric. Notice that 
only the first three terms in the action above contain linear perturbations around the chosen FRW 
solution, all the others are explicitly quadratic or higher. Therefore the coefficients c{t) and h.{t) 
will be fixed by the requirement of having a given FRW evolution H{t), i.e. requiring that tadpole 
terms cancel around this solution. Before fixing these coefficients, it is important to realize that 
this simplification is not trivial. One would expect that there are an infinite number of operators 
which give a contribution at first order around the background solution. However one can write 
the action as a polynomial of linear terms like SK^^i, and g^^ + 1, so that it is evident whether an 
operator starts at linear, quadratic or higher order. All the linear terms besides the ones in eq. ([5]) 
will contain derivatives and they can be integrated by parts to give a combination of the three linear 
terms we considered plus covariant terms of higher order. This construction is explicitly carried out 
in appendix iBl We conclude that the unperturbed history fixes c{t) and A{t), while the difference 
among different models will he encoded into higher order terms. 

We can now fix the linear terms imposing that a given FRW evolution is a solution. As we 
discussed, the terms proportional to c and A are the only ones that give a stress energy tensor 

which does not vanish at zeroth order in the perturbations and therefore contributes to the right 
hand side of the Einstein equations. During inflation we are mostly interested in a flat FRW Universe 
(see Appendix [B] for the general case) 

ds'^ = -dt^ ^a^[f)dx^ (7) 

so that Friedmann equations are given by 

^'^^'^^ = -^N*)-A(*^ (9) 

Solving for c and A we can rewrite the action ([5]) as 

S = id^x \/^[^M|,fi + M|,«9»» - M2,{3S2 + Jj) + iMjlt)"!;,"" + 1)^ + \^Kh[{)\g'« + 1)" + 
-^(,» + l)„^<.„ - Mt,K^^^ _ ^^A-„^A-„ + ...] . (10) 

As we said all the coefficients of the operators in the action above may have a generic time depen- 
dence. However we are interested in solutions where K and K do not vary significantly in one Hubble 
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time. Therefore it is natural to assume that the same holds for all the other operators. With this 
assumption the Lagrangian is approximately time translation invariant [§. Therefore the time depen- 
dence generated by loop effects will be suppressed by a small breaking parameter H. This assumption 
is particularly convenient since the rapid time dependence of the coefficients can win against the fric- 
tion created by the exponential expansion, so that inflation may cease to be a dynamical attractor, 
which is necessary to solve the homogeneity problem of standard FRW cosmology. 

It is important to stress that this approach does describe the most generic Lagrangian not only 
for the scalar mode, but also for gravity. High energy effects will be encoded for example in operators 
containing the perturbations in the Riemann tensor SR^ypa. As these corrections are of higher order 
in derivatives, we will not explicitly talk about them below. 

Let us give some examples of how to write simple models of inflation in this language. A model 
with minimal kinetic term and a slow-roll potential V{4>) can be written in unitary gauge as 



- V{<t>) 



(fx 



(11) 



As the Friedmann equations give 4>o{tf = -2Mj,H and V{(p{t)) = M^y{3H^ + H) we see that the 
action is of the form (jlOp with all but the first three terms set to zero. Clearly this cannot be true 
exactly as all the other terms will be generated by loop corrections: they encode all the possible 
effects of high energy physics on this simple slow-roll model of inflation. 

A more general case includes all the possible Lagrangians with at most one derivative acting on 
each (p: L = P{X, 0), with X 



g^'^df^(f)d,^(j). Around an unperturbed solution </>o(t) we have 
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d^xV^P(<Ao(t)V°,<A(i)) 



which is clearly of the form above with M^{t) = (j)o{t)'^"'d"'P/dX'^ evaluated at (j)o{t). Terms con- 
taining the extrinsic curvature contain more than one derivative acting on a single scalar and will 
be crucial in the limit of exact de Sitter, H ^ 0. 



3 Action for the Goldstone Boson 

The unitary gauge Lagrangian describes three degrees of freedom: the two graviton helicities and 
a scalar mode. This mode will become explicit after one performs a broken time diffeomorphism 
(Stiickelberg trick) as the Goldstone boson which non-linearly realizes this symmetry. In analogy 
with the equivalence theorem for the longitudinal components of a massive gauge boson [9], we expect 

^The limit in which the time shift is an exact symmetry must be taken with care because H ^ 0. This 
implies that the spatial kinetic term for the Goldstone vanishes, as we will see in the discussion of Ghost 
Inflation. 

•^Notice that this symmetry has nothing to do with the breaking of time diffeomorphisms. To see how this 
symmetry appears in the (j) language notice that, after a proper field redefinition, one can always assume that 
cj) = const. With this choice, invariance under time translation in the unitary gauge Lagrangian is implied by 
the shift symmetry cj) (j) + const. This symmetry and the time translation symmetry of the (j) Lagrangian 
are broken down to the diagonal subgroup by the background. This residual symmetry is the time shift in 
the unitary gauge Lagrangian. 
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that the physics of the Goldstone decouples from the two graviton hehcities at short distance, when 
the mixing can be neglected. Let us review briefly what happens in a non-Abelian gauge theory 
before applying the same method in our case. 



The unitary gauge action for a non-Abelian gauge group is 



S 



where = A^T'^. Under a gauge transformation we have 

A, 

The action therefore becomes 



UAM^ + -UdJJ^ = -UDM^ . 
9 9 



1 



1 m 
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2 TtD^U^D^U . 



(13) 



(14) 



(15) 



The gauge invariance can be "restored" writing U = exp[iT""rr'^{t,x)], where vr" are scalars (the 
Goldstones) which transform non-linearly under a gauge transformation A as 



(16) 



Going to canonical normalization tTc = m/g ■ ir, we see that the Goldstone boson self-interactions 
become strongly coupled at the scale 4:7rm/g, which is parametrically higher than the mass of the 
gauge bosons. The advantage of reintroducing the Goldstones is that for energies E ^ m the mixing 
between them and the transverse components of the gauge field becomes irrelevant, so that the two 
sectors decouple. Mixing terms in eq. ()14p are in fact of the form 



m 



(17) 



which are irrelevant with respect to the canonical kinetic term {diTc)^ for E ^ m. In the window 
m <^ E <^ AiTm/g the physics of the Goldstone vr is weakly coupled and it can be studied neglecting 
the mixing with transverse components. 

Let us follow the same steps for our case of broken time diffeomorphisms. Let us concentrate for 
instance on the two operators: 



j <fx^g[A{t) + B{t)g''\x)] . 

Under a broken time diff. t ^ t = t + C^{x), x ^ x = x, g^^ transforms as: 

00/ N ~O0/~/ N\ dx^{x) dx^{x) 
g ix) g {x(x)) = — g^ (x) . 

The action written in terms of the transformed fields is given by: 

dx^ dx^ 



(18) 



d!^x \/ —g{x{x)) 



dx 




dx 





A{t)+B{t) 



dxt^ dx 



^r^ixix)) 



(19) 



(20) 
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Changing integration variables to x, we get: 



A{t-eixix)))+Bit-eix{x))) 



dit - ^^{x{x))) dit - ^%xix))) 



dxi' 



dx'' 



(21) 



The procedure to reintroduce the Goldstone is now similar to the gauge theory case. Whenever 
appears in the action above, we make the substitution 



fix{x)) -tt(x) . 



This gives, dropping the tildes for simplicity: 



/ 



d^x y/-g{x) 



A{t + 7r{x)) + B{t + TT{x)) 



d{t + 7r{x)) d{t + Tr{x)) 



dxi^ 



dx'^ 



g^'ix) 



(22) 



(23) 



One can check that the action above is invariant under diffs at all orders (and not only for infinitesimal 
transformations) upon assigning to vr the transformation rule 



tt{x) —>■ 7r(x(x)) = 7r(x) — £,^{x) . 

With this definition tt transforms as a scalar field plus an additional shift under time diffs. 
Applying this procedure to the unitary gauge action ([10]) we obtain 



(24) 



S 



(fx 
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M^Am\t + -k) + H{t + -k)] + 



(25) 



Ml^H{t + vr) ((1 + ^)2g00 + 2(1 + iT)diTTg''' + g'^diTrd^Tr) + 



M2{t 



2! 

M3(t + 7r)4 



3! 



((1 + 7^)2500 + 2(1 + 7^)ai7rr7°* + g'^ diirdj^ + l) + 
((1 + 7^)2500 + 2(1 + 7^)9,7150* + g'^diT^djTi + + ... 



where for the moment we have neglected for simplicity terms that involve the extrinsic curvature. 

This action is rather complicated, and at this point it is not clear what is the advantage of 
reintroducing the Goldstone vr from the unitary gauge Lagrangian. In analogy with the gauge 
theory case, the simplification occurs because, at sufficiently short distances, the physics of the 
Goldstone can be studied neglecting metric fluctuations. As for the gauge theory case, the regime 
for which this is possible can be estimated just looking at the mixing terms in the Lagrangian above. 
In eq. (125p we see in fact that quadratic terms which mix vr and g^y contain fewer derivatives than 
the kinetic term of vr so that they can be neglected above some high energy scale. In general the 
answer will depend on which operators are present. Let us start with the simplest case in which only 
the tadpole terms are relevant (M2 = M3 = ... = 0). This corresponds to the standard slow-roll 
inflation case. The leading mixing with gravity will come from a term of the form 



M§^Hi:6g 



00 



(26) 



After canonical normalization (vTc ~ Mp\H^^'^iT, 5g^ ~ M-pidg^'^), we see that the mixing terms 
can be neglected for energies above -Emix ~ e^/'^H, where e is the usual slow-roll parameter e = 
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-H/H'^. Another case which will be of interest is when the operator M2 gets large. In this case we 
have mixing terms of the form 

~ M|^5ff°° (27) 

which, upon canonical normalization (notice that now tTc ~ M|7r), becomes negligible at energies 
larger than ^mix ~ M|/Mpi 0. 

In the regime E ^ -Emix the action dramatically simplifies to 



.(28) 



Given an inflationary model, one is interested in computing predictions for present cosmological 
observations. From this point of view, it seems that the decoupling limit (j28|) is completely irrelevant 
for these extremely infrared scales. However, as for standard single field slow-roll inflation, one can 
prove that there exists a quantity, the usual C, variable, which is constant out of the horizon at any 
order in perturbation theory [10^ [TT] (see Appendix D of [12j for a generalization including terms 
with higher spatial derivatives). The intuitive reason for the existence of a conserved quantity is 
that after exiting the horizon different regions evolve exactly in the same way. The only difference is 
how much one has expanded with respect to another and it is this difference that remains constant. 

Therefore the problem is reduced to calculating correlation functions just after horizon crossing. 
We are therefore interested in studying our Lagrangian with an IR energy cutoff of order H. If the 
decoupling scale -Emix is smaller than H, the Lagrangian for vr (128^ will give the correct predictions 
up to terms suppressed by E^ix/H. 

As we discussed, we are assuming that the time dependence of the coefficients in the unitary 
gauge Lagrangian is slow compared to the Hubble time, that is, suppressed by some generalized slow 
roll parameters. This implies that the additional vr terms coming from the Taylor expansion of the 
coefficients are small. In particular, the relevant operators, i.e. the ones which dominate moving 
towards the infrared, like the cubic term, are unimportant at the scale H and have therefore been 
neglected in the Lagrangian ([28|) . 

In conclusion, with the Lagrangian ()28p one is able to compute all the observables which are not 
dominated by the mixing with gravity, like for example the non-Gaussianities in standard slow-roll 
inflation [131 114j . Notice however that the tilt of the spectrum can be calculated, at leading order, 
with the Lagrangian (f28ll . As we will see later, its value can in fact be deduced simply by the 
power spectrum at horizon crossing computed neglecting the mixing terms. It is important to stress 
that our approach does not lose its validity when the mixing with gravity is important so that the 
Goldstone action is not sufficient for predictions. The action (jlOp contains all the information about 
the model and can be used to calculate all predictions even when the mixing with gravity is large. 



"^In the theories we are studying Lorentz symmetry is spontaneously broken, so one should define a separate 
regime of energies and momenta for which the mixing can be neglected. For cosmological perturbations, we 
will be only interested in the energy range. 
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4 The various limits of single field infiation 



4.1 Slow-roll inflation and high energy corrections 

The simplest example of the general Lagrangian (jlOp is obtained by keeping only the first three 
terms, which are fixed once we know the background Hubble parameter H{t), and setting to zero 
all the other operators of higher order: M2 = M3 = Mi = M2 . . . = 0. In the (j) language, this 
corresponds to standard slow-roll inflation, with no higher order terms. In this discussed 
in the last section, predictions at the scale H can be made neglecting the mixing with gravity and 
concentrating on the Goldstone Lagrangian (|28p . One is interested in calculating, soon after horizon 
crossing, the conserved quantity C- This is defined, at linear order, by choosing the gauge vr = 
(unitary gauge in our language) and the spatial part of the metric to be 

gij = a^t) [(1 + 2C(t, x))5ij + -fij] (29) 

where 7 is transverse and traceless and it describes the two tensor degrees of freedom. The relation 
between vr and C is very simple. As we are neglecting the mixing with gravity, the metric is unper- 
turbed in the vr language; to set vr = one has to perform a time diffeomorphism t ^ t — 7r(t, x) 
which gives a spatial metric of the form (j29[) with 

C{t,x) = -HTT{t,x) . (30) 

For each mode k, one is only interested in the dynamics around horizon crossing io{k) = k/a ^ H . 
During this period the background can be approximated as de Sitter up to slow-roll corrections. 
Therefore, the 2-point function of the canonically normalized scalar tTc is given by the de Sitter 
result 



(7r,(A:i)7r,(A:2)) = (27r)^,5(A;i + A^a)^ 



tt2 

(31) 



where here and below * means the value of a quantity at horizon crossing. This implies that the 
2-point function of C is given by 

rnxm = (2,mk. + = i^-'M^^ + &)^^ . (32) 

As the variable ( is constant outside the horizon, this equation is exact for all k up to slow-roll 
corrections. In particular it allows us to calculate the tilt of the spectrum at leading order in 
slow-roll 

1 d 1 d Hf 

rig — I = — r log —. — = -r— ^— log —. — = 4-— 77 — . (33) 

dlogk H,dt, H,H, ^ ' 

Notice however that not all observables can be calculated from the vr Lagrangian (|28p : this 
happens when the leading result comes from the mixing with gravity or is of higher order in the 
slow-roll expansion. For example, as the first two terms of eq. (|28p do not contain self-interactions of 
TT, the 3-point function {({ki)({k2)C{k3)) would be zero. One is therefore forced to look at subleading 
corrections, taking into account the mixing with gravity in eq. (j25p . 
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Obviously our choice of setting to zero all the higher order terms cannot be exactly true. At the 
very least they will be radiatively generated even if we put them to zero at tree level. The theory is 
non-renormalizable and all interactions will be generated with divergent coefficients at sufficiently 
high order in the perturbative expansion. As additional terms are generated by graviton loops, 
they may be very small. For example it is straightforward to check that starting from the unitary 
gauge interaction Mp^Hg^^ a term of the form [g^^ + 1)^ will be generated with a logarithmically 
divergent coefficient ~ iJ^logA. This implies that one should assume M| > (jl). This 
lower limit is however very small. For example the dispersion relation of vr will be changed by 
the additional contribution to the time kinetic term: this implies, as we will discuss thoroughly 
below, that the speed of vr excitations deviates slightly from the speed of light, by a relative amount 
1 — Cs ~ M2 / {\H\Mpy) ~ \H\/Mp^. Using the normalization of the scalar spectrum eq. (j32p . we see 
that the deviation from the speed of light is > • 10^^^. A not very interesting lower limit. 

The size of the additional operators will be much larger if additional physics enters below the 
Planck scale. In general our approach gives the correct parametrization of all possible effects of 
new physics. As usual in an effective field theory approach, the details of the UV completion of 
the model are encoded in the higher dimension operators. This is very similar to what happens 
in physics beyond the Standard Model. At low energy the possible effects of new physics are 
encoded in a series of higher dimensional operators compatible with the symmetries |16j . The 
detailed experimental study of the Standard model allows us to put severe limits on the size of these 
higher dimensional operators. The same can be done in our case, although the set of conceivable 
observations is unfortunately much more limited. One example of a possible experimental limit on 
higher dimension operators is the consistency relation for the gravitational wave tilt. As is well 
known, the gravity wave spectrum from the Einstein-Hilbert action is given by 

{Y{ki)Y'{k2)) = i27rf6{h + ^2)^^<5s.' (34) 

where 7^ denotes the two possible polarizations of the gravity wave. The ratio between this contri- 
bution and the scalar one (j32p is given by e^=. The gravitational wave tilt, Ug = — 2e^,, is thus fixed 
once the ratio between tensor and scalar modes is known. 

This prediction is valid if one assumes M2 = 0, i.e. Cg = 1. As we will see in fact, the 
scalar spectrum goes as cj^, while predictions for gravitational waves are not changed by M2. The 
experimental verification of the consistency relation, even with large errors, would tell us that Cg 
cannot deviate substantially from 1 which implies 

M| < M^,\H\ . (35) 

Notice that the higher dimension operators will not only influence scalar fluctuations, but also the 
tensor modes, although these corrections are arguably much harder to test. For example the unitary 
gauge operator —M'i{t)/2-5K^y5K^^, whose relevance for scalar fluctuations will be discussed later 
on, contains terms of the form gj-. This will change the gravity wave dispersion relation. It is in 
fact straightforward to obtain the action for the tensor modes 7jj in the presence of this operator. 

^The explicit calculation of logarithmic divergences in a theory of a massless scalar coupled to gravity has 
been carried out a long time ago in [15j . 
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One gets 



8 



-g 



1 



Ml, 



(36) 



Therefore the spectrum of gravity waves (I34p will get corrections of order M|/Mpj. This correction 
is small unless we push Mg up to the Planck scale. It is easy to realize that operators of the form 
{g^^ + !)"■ do not influence tensor modes as they do not affect the transverse-traceless components 
of the metric. 

Other examples of experimental limits on various operators will be discussed in the following 
sections. 



4.2 Small speed of sound and large non-Gaussianities 

The Goldstone action (1281) shows that the spatial kinetic term {diir)^ is completely fixed by the 
background evolution to be Mp,H{diTT)^. In particular only for H < 0, it has the "healthy" negative 
sign. This is an example of the well studied relationship between violation of the null energy 
condition, which in a FRW Universe is equivalent to iiT < 0, and the presence of instabilities in the 
system [HI [8]. Notice however that the wrong sign of the operator (5j7r)^ is not enough to conclude 
that the system is pathological: higher order terms like 5K^^ ^ may become important in particular 
regimes, as we will discuss thoroughly below. Reference [1] studies examples in which H > can be 
obtained without pathologies. 

The coefficient of the time kinetic term vr^ is, on the other hand, not completely fixed by the 
background evolution, as it receives a contribution also from the quadratic operator {g^^ + 1)^. In 
eq. (|28ll we have 

' -M^iij + 2M|) 7i-2 . (37) 



To avoid instabilities we must have —Mp^H + 2M2 > . As time and spatial kinetic terms have 
different coefficients, vr waves will have a "speed of sound" Cs 7^ 1. This is expected as the background 
spontaneously breaks Lorentz invariance, so that = 1 is not protected by any symmetry. As we 
discussed in the last section, deviation from = 1 will be induced at the very least by graviton 
loops [§. The speed of sound is given by 

= 1 ^ . 38 

This implies that in order to avoid superluminal propagation we must have M| > (assuming H < 
0). Superluminal propagation would imply that the theory has no Lorentz invariant UV completion 
|18j . In the following we will concentrate on the case Cg < 1, see [19] for a phenomenological 
discussion of models with c, > 1. 



^If we neglect the coupling with gravity and the time dependence of the operators in the unitary gauge 
Lagrangian (so that n ^ n + const is a symmetry), Cg = 1 can be protected by a symmetry 9p7r S^tt + v^, 
where i;^ is a constant vector. Under this symmetry the Lorentz invariant kinetic term of tt changes by a 
total derivative, while the operator proportional to M| in eq. ([28)1 is clearly not invariant, so that — 1. 
Notice that the theory is not free as we are allowed to write interactions with more derivatives acting on tt. 
This symmetry appears in the study of the brane bending mode of the DGP model jl8| . 
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Using the equation above for the Goldstone action can be written at cubic order as 



.(39) 



From the discussion in section ([3]) we know that the mixing with gravity can be neglected at 
energies E ^ -Emix — AffZ-^Pi- This imphes that predictions for cosmological observables, which 
are done at energies of order H, are captured at leading order by the Goldstone action if 
H » M|/Mpi, or equivalently for e/c^ ^ 1. If this is not the case one is not assured that the 
Goldstone action contains the leading effects. 

The calculation of the 2-point function follows closely the case c<j = 1 if we use a rescaled 



momentum k = Cgk and take into account the additional factor c<, ^ in front of the time kinetic term. 
We obtain 



(C(^i)C(4)) = {27rfS{h + h)— ■ TTT^ A = {27Tf5{h + ^2)— • y^TFl^ • (40) 



The variation with time of the speed of sound introduces an additional contribution to the tilt 

n, = ^log^^ = ±±logJ^= - ^ . (41) 



The result agrees with the one found in |20j . 

From the action (j39p we clearly see that the same operator giving a reduced speed of sound 
induces cubic couplings of the Goldstones of the form 7T(V7r)^ and vr'^. The non- linear realization of 
time diffeomorphisms forces a relation between a reduced speed of sound and an enhanced level of 
the 3-point function correlator, i.e. non-Gaussianities. This relationship was stressed in the explicit 
calculation of the 3-point function in [21] . 

To estimate the size of non-Gaussianities, one has to compare the non-linear corrections with 
the quadratic terms around freezing, uj ^ H. In the limit Cs <C 1, the operator 7T(V7r)^ gives the 
leading contribution, as the quadratic action shows that a mode freezes with k ~ H/cs., so that 
spatial derivatives are enhanced with respect to time derivatives. The level of non-Gaussianity will 
thus be given by the ratio: 



-7i-(V7r)2 



^vr-^C, (42) 



where in the last step we have used the linear relationship between vr and C, eq. ()30|) . Taking 
^ ~ 10~^ we have an estimate of the size of the non-linear correction 0. Usually the magnitude of 



^The size of the non-linear corrections depend on the specific value of Even if the typical value of C, is 
small, one may be interested in very large (and therefore very unlikely) fluctuations, for example to study the 
production of primordial black holes. For sufliciently large values of C, C > Cg, non-linear corrections become 
of order 1 and the perturbative expansion breaks down. Therefore, predictions which depend on very large 
values of ^ may lie out of the regime of validity of the effective field theory. 
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non-Gaussianities is given in terms of the parameters /nlj which are parametrically of the form: 
'^7i-(V7r)2/^2 ~ /nlC- The leading contribution will thus give 

/NL!'7r{V7r)2 ~ 3 ' (43) 

The superscript "equil." refers to the momentum dependence of the 3-point function, which in these 
models is of the so called equilateral form [22] . This is physically clear in the Goldstone language as 
the relevant vr interactions contain derivatives, so that they die out quickly out of the horizon; the 
correlation is only among modes with comparable wavelength. 

In the Goldstone Lagrangian (I39p there is an additional independent operator, —^M^tt^, con- 
tributing to the 3-point function, coming from the unitary gauge operator (g^^ + 1)^. We thus have 
two contributions of the form tt^ which give 

j^equil. , 4 M| 

^NL-^^'-3l^|^- 

The size of the operator —^M^n^ is not constrained by the non-linear realization of time diffeomor- 
phisms: it is a free parameter. In DBI inflation [2] we have M| ~ Mpj|if|c~^, so that its contribution 
to non-Gaussianities is of the same order as the one of eq. (j43p . The same approximate size of the 
-M| is obtained if we assume that both the unitary gauge operators M|(g^^ + 1)^ and M^{g^^ + 1)^ 
become strongly coupled at the same energy scale. 

It is interesting to look at the experimental limits on non-Gaussianities as a constraint on the 
size of the unitary gauge operator {g^^ + 1)^ and therefore on the speed of sound. The explicit 
calculation [21j gives the contribution of the operator niVn)'^ to the experimentally constrained 
parameter /i^*?"'''; at leading in order in we have[f| 



/r'=ii. (45) 



The experimentally allowed window 

- 256 < /X''' < 332 at 95% C.L. (46) 

translates into the constraint 



Cs > 0.028 at 95% C.L. (47) 

Notice however that, although in principle the operators 7r(V7r)2 and vr^ give a different momentum 
dependence of the 3-point function, this difference is not experimentally appreciable at present, so 
that the constraint (j46p is on the joint effect of the two operators. The constraint on the speed of 
sound will hold barring a cancellation between the two operators. In the case of DBI inflation for 
example the effect of the operator M^{g^^ + 1)^ is sizeable as we discussed. However there is no 
cancellation and the constraint on the speed of sound is only slightly changed to 

DBI : Cs > 0.031 at 95% C.L. (48) 
®This is obtained setting P^xxx = in the notation of [?!] . 
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Although we concentrated so far on the Goldstone Lagrangian, it is important to stress that this 
general approach is useful also when one is interested in taking into account the full mixing with 
gravity. For example, going back to the unitary gauge Lagrangian (jlOp . we can easily see how many 
coefficients will be relevant in calculating the 3-point function. At leading order in slow-roll and in 
derivatives there are 2 coefficients as we discussed: M2 and M3. At first order in slow-roll, there 
will be 4 new parameters describing the slow variation of the coefficients: the conventional e and 
•q slow-roll parameters and two additional ones for the coefficients of the operators {g^^ -|- 1)^ and 
{g^^ + 1)^- This in fact is what one finds in the explicit calculation |21j |^. 

All the discussion can be straightforwardly extended to the 4-point function (and higher order 
correlators). In the Goldstone Lagrangian we have 3 operators contributing to the 4-point function 
(again at leading order in slow-roll and derivatives): {g^^ + 1)^, which is fixed by the speed of sound 

{9^^ + 1)^ and {g^^ + f)^. Let us estimate the effect of the operator which is fixed by the speed 
of sound. As we did for the 3-point function, it is easy to see that the effect will be dominated by 
the operator (Vvr)^ and that the level of non-Gaussianity induced by it can be estimated as 

W JS-)\,ff%. 1^., (49) 



C2 H^TT^ Ci C\ 

This matches with the explicit calculation done in 



4.2.1 Cutoff and Naturalness 

As discussed, for < 1 the Goldstone action contains non-renormalizable interactions. Therefore 
the self-interactions among the Goldstones will become strongly coupled at a certain energy scale, 
which sets the cutoff of our theory. This cutoff can be estimated looking at tree level partial wave 
unitarity, i.e. finding the maximum energy at which the tree level scattering of vrs is unitary. The 
calculation is straightforward, the only complication coming from the non-relativistic dispersion 
relation. The cutoff scale A turns out to be 

7 5 

A^ ~ 16^2^1—^^ ~ X^T^Ml^T^ ■ (50) 
(1 - ciY l-cj 

The same result can be obtained looking at the energy scale where loop corrections to the vrvr 
scattering amplitude become relevant. As expected the theory becomes more and more strongly 
coupled for small Cg, so that the cutoff scale decreases. On the other hand, for — > 1 the cutoff 
becomes higher and higher. This makes sense as there are no non-renormalizable interactions in this 
limit and the cutoff can be extended up to the Planck scale. This cutoff scale is obtained just looking 
at the unitary gauge operator {g^^ + 1)^; depending on their size the other independent operators 
may give an even lower energy cutoff. Notice that the scale A indicates the maximum energy at 
which our theory is weakly coupled and make sense; below this scale new physics must come into 
the game. However new physics can appear even much below A. 

If we are interested in using our Lagrangian for making predictions for cosmological correlation 
functions, then we need to use it at a scale of order the Hubble parameter H during inflation. We 

^The explicit calculation shows that one of the coefficients does not give rise to an independent momentum 
dependence of the 3-point function, so that it cannot be disentangled from the other parameters. 
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therefore need that this energy scale is below the cutoff, H <^ A. Using the explicit expression for 
the cutoff ([50]) in the case ^ 1 one gets 

H"^ < M^^\H\4 (51) 

which can be rewritten using the spectrum normalization (j40p as an inequality for the speed of sound 

Cs > P^^"^ ~ 0.003 . (52) 

A theory with a lower speed of sound is strongly coupled &t E H. Not surprisingly this value 
of the speed of sound also corresponds to the value at which non-Gaussianity are of order one: the 
theory is strongly coupled at the energy scale H relevant for cosmological predictions. 

Let us comment on the naturalness of the theory. One may wonder whether the limit of small 
Cs is natural or instead loop corrections will induce a larger value. The Goldstone self- interactions, 
7r(V7r)^ and (Vvr)^ for example, will induce a radiative contribution to (Vvr)^. It is easy to estimate 
that these contributions are of order c7^A^/(167r^M|), where A is the UV cutoff, i.e. the energy scale 
at which new physics enters in the game. We can see that it is impossible to have large radiative 
contribution; even if we take A at the unitarity limit (j50p . the effect is of the same order as the 
tree level value. This makes sense as the unitarity cutoff is indeed the energy scale at which loop 
corrections become of order one. 

We would like also to notice that the action (128p is natural from an effective field theory point of 
view |28j . The relevant operators are in fact protected from large renormalizations if we assume an 
approximate shift symmetry of vr. In this case the coefficients of the relevant operators are sufficiently 
small and they will never become important for observations as cosmological correlation functions 
probe the theory at a fixed energy scale of order H: we never go to lower energy. Clearly here we 
are only looking at the period of inflation, where an approximate shift symmetry is enough to make 
the theory technically natural; providing a graceful exit from inflation and an efficient reheating are 
additional requirements for a working model which are not discussed in our formalism. 



4.3 De-Sitter Limit and the Ghost Condensate 

In the previous section we saw that the limit Cs — > is pathological as the theory becomes more and 
more strongly coupled. However we have neglected in our discussion the higher derivative operators 
in the unitary gauge Lagrangian (fTO]) 



-9 



M2{tf 



(53) 



These operators give rise in the Goldstone action to a spatial kinetic term of the form 



1 



(Sfvr) 



(54) 



where = M| + M|. This spatial kinetic term will make the Goldstone propagate even in the 
limit Cs 0. It is therefore interesting to consider our general Lagrangian in the limit H = 0, when 
the gravitational background is exactly de Sitter space which implies Cs = 0. As H is now time 
independent, it is possible to impose an additional symmetry to the theory: the time independence 
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of all the coefficients in the unitary gauge Lagrangian. Looking back at the procedure (j25p to 
reintroduce the Goldstone vr, we realize that this symmetry forbids any dependence on vr without 
derivatives. The Goldstone action is thus invariant under shift of vr 



vr(x, t) vr(x, t) + const. (55) 

This is the limit of Ghost Condensation [25], where the Goldstone has a non-relativistic dispersion 
relation lo (x k"^. More generally one can consider intermediate situations where both the spatial 
kinetic term Cg(Vvr)^ and the higher derivative one (V^vr)^ are present. The predictions of the theory 
will change significantly depending on which term dominates at the energy of freezing oj ^ H [71 127j . 

As with the previous models, one must find the energy regime for which the mixing of the 
Goldstone with gravity can be neglected. One simple way to estimate this range is to look at the 
5K operators which contain terms like 

5Kij D {didj-K + dtgoj) . (56) 

Going to canonical normalization this shows that the mixing with gravity can be neglected for 
k > M|/Mpi. As the dispersion relation of the Goldstone is of the form o;^ = M^/M| • /c^, we see 
that the energy -Emix under which the mixing is relevant is -Emix — MM2/Mpy [25]. Notice that this 
scale has nothing to do with the the curvature of the background. This is a quite remarkable feature 
of this example, as usually the mixing with gravity is related to the background stress energy tensor 
and therefore to the curvature of spacetime: the more a system curves space, the more it mixes 
with gravity. In this case on the other hand, the mixing will be relevant even on a flat Minkowski 
background. This is what one calls a proper modification of gravity: gravity, for example the 
Newtonian potential generated by a source, is modified at scales much smaller than the curvature. 
This model of modification of gravity and its rich phenomenology has been studied in [25j 0. 

As we are interested in inflation, we concentrate on the opposite limit H S> -Emix) when the mixing 
can be neglected and one can focus on the vr Lagrangian. Let us briefly describe the main features of 
Ghost Inflation, referring for details to [251 [71 [27]. where the theory is studied with an approach very 
close to the one presented in this paper. Most of the interesting features can be understood looking 
at the scaling with energy of the various operators. Given the non-relativistic dispersion relation, 
uj (X k'^, the way an operator scales with energy does not coincide with its mass dimension as in 
the Lorentz invariant case. A rescaling of the energy by a factor s, E ^ sE, (equivalent to a time 
rescaling t — > s~^t), must go together with a momentum transformation k s^^'^k (x s~^^'^x on 
the spatial coordinates). As the quadratic action for the Goldstones is of the form 



2M,V-^l(a^) 



(57) 



we must assign to vr the scaling dimension 1/4 

vr^s^/V. (58) 



^°Also in the case of models with a reduced speed of sound, the scale of mixing with gravity can become 
parametrically smaller than the horizon; it is enough to have e/c^ ^1. In this case the model can be 
considered a way of modifying gravity. Notice however that one can not take the limit H ^ without 
considering the spatial higher derivative terms: the scalar mode would not propagate otherwise. 
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to keep the quadratic action invariant. With this rule it is easy to check that ah the ahowed 
Goldstone operators, besides the kinetic term (j57p . are irrelevant, i.e. they have positive scaling 
dimension and they become less and less relevant going down in energy. This shows that the theory 
makes sense as an effective field theory. In particular the higher derivative time kinetic operator 
71"^, which would naively seem as important as (V^vr)^ and would describe the presence of a ghost 
in the theory, has dimension 2 and it can be neglected at low energies. If one assumes that there is 
a single scale M in the problem, M ~ M2 — M, this will also set the energy cutoff of the effective 
field theory description. 

The scaling dimension of vr also allows us to estimate the spectrum of perturbations produced in 
Ghost Inflation. The dimension of vr tells us how the amplitude of quantum fluctuations changes with 
energy. At the scale of the cutoff A ~ M, the quantum fluctuations of the canonically normalized 
Goldstone field tTc — M^vr are of the order of the cutoff S7rc{M) ~ M. Going down in energy we can 
estimate the quantum fluctuations at freezing E ^ H. In the standard case the scalar would have 
dimension 1 and its fluctuations at freezing would be of order H] in this case on the other hand we 
have 

/ tr \ 1/4 

5n,{H) ~ 5^,{M) i — \ ~ {HmY^ . (59) 

Quantum fluctuations at the scale H are much enhanced with respect to a scalar with a Lorentz 
invariant dispersion relation. The spectrum of C will thus be given by [7] 

(C(^i)C(^2)) ~ {27,f6{h + k2)^,{HM')y'^^ . (60) 

The correct normalization of the spectrum requires 

The non-Gaussianity will be dominated by the operator with the lowest dimension. It is straight- 
forward to see that the operator M27r(V7r)^ coming from the unitary gauge operator M2 (g'^'^ + l)^ has 
dimension 1/4 and it is the least irrelevant operator. At the cutoff scale M the theory is strongly 
coupled. As the cubic operator has dimension 1/4, at energies of order H it will give a level of 
non-Gaussianity of order (H/M)^/^, which is parametrically of order P^^^ . The same result can be 
obtained with the approach used in the last section, i.e. comparing the interaction term with the 
free action at freezing 

^ M-H.jHM).^ = = ~ f-V' (61) 
C2 M4i/V2 \m) • 

The level of non-Gaussianity is extremely high compared to standard slow-roll as a consequence of 
the very low dimension of the most relevant operators. The explicit calculation [7] gives an effect 
which is somewhat smaller than the naive estimate and comparable to the existing experimental 
bound [23]. 

In our discussion we have neglected so far the unitary gauge operator 

d'x (-^^{9'' + l)dK^,) . (62) 



This operator is odd under time reversal, so that it is consistent to set it to zero. If this term is 
present, there is a second operator with dimension 1/4 in the Goldstone Lagrangian, of the form 
V'^tt('Vtt)'^ . Its contribution to the 3-point function would be comparable with 7r(V7r)2. The unitary 
gauge operator (f62l) also contributes to the vr quadratic Lagrangian as we are now going to discuss. 
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4.3.1 De-Sitter Limit without the Ghost Condensate 



In this section we want to study the effect of the operator (|62|) on the quadratic vr action. We will 
see that, if the coefficient of this operator is sufficiently large, we obtain a new de Sitter limit, where 
the dispersion relation at freezing is of the form uj^ cck'^ , instead of the Ghost Condensate behavior 
oc k*. 

For simplicity we can take Mi to be time independent. Reintroducing the Goldstone we get a 
3-derivative term of the form — Mj'vrV'^Tr/a^ (0). This would be a total time derivative without the 
time dependence of the scale factor a{t) and of the metric determinant. Integrating by parts we get 
a standard 2-derivative spatial kinetic term 



/ 



rf.,^^*!^/'^., , (63) 



a 



In the exact de Sitter limit, H = 0, and taking M2 ~ Mi ~ M, this operator gives a dispersion 
relation of the form o;^ = c^/c^, with a small speed of sound 

= — <C 1 . (64) 

This will hold only if the higher derivative operators 5K^ ^ and SK'^^dK'^ are subdominant. If we 
assume that they are characterized by the same mass scale, M2 ~ M3 ~ M, the dispersion relation 
will get two contributions 

The two spatial kinetic term are comparable at freezing ui ^ H. On the other hand, if the k^ 
contribution is somewhat suppressed, it becomes irrelevant at freezing and therefore for inflationary 
predictions. In this limit we have a new kind of Ghost Inflation with an exactly de Sitter background, 
but with a cj^ oc fc^ dispersion relationship at freezing. 

Following what we did for the other models it is straightforward to obtain the spectrum normal- 
ization and an estimate of the 3-point function non-Gaussianity. 

(C(fci)C(^2)) - {27Tf6{ki + k)T7I^T^- {27rf6{h +h) tj) 71- (66) 



M'^cjkf ^ ' ' ^'\M) kf 



\Hnr.(^Y\ (67) 



£2 ^ c2"" Km J 



-'^-'^The operator gives also a contribution to tt^ proportional to H. We will assume that this is small compared 
to MjTT^. In Minkowski space the operator we are studying can be forbidden by a </) ^ ^(j) symmetry, which 
is equivalent to time reversal in unitary gauge [25j . In a de Sitter background this symmetry is broken by the 
metric, so that this operator cannot be set to zero. 

-"^^In this model the mixing with gravity is rather different from the previous cases. The reason is that a 
time derivative is integrated by parts to get to eq. (|63p , so that the Goldstone terms contain the same number 
of derivatives as the terms describing the mixing with gravity. This implies that the mixing does not become 
less and less relevant going to high energy. On the other hand one can choose the model parameters in such 
a way that the mixing is always irrelevant. See |23j for the explicit calculations. 
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A comparable contribution will come from the Goldstone operator (V7r)^V^7r. Not surprisingly the 
estimates above are the same as the ones we obtained in the Ghost Condensate case eqs (|60p and 
(j6ip . As we discussed in fact, taking all the operators at the same scale one gets a comparable 
contribution at freezing from the k'^ and k^ spatial kinetic terms. We thus expect similar predictions 
when we assume that only one of the two contributions is present. 

Now that we have found two different de Sitter limits, one dominated at freezing by {g^^ + l)5K^^^ 
and the other by dRf"^^ and SK^'^dK"^, one may wonder if there are other possibilities. One 
could imagine that both these spatial kinetic terms are suppressed for some reason and the leading 
operators come at higher order. In this case one would end up with a dispersion relation of the form 

a;2 ^ A;2" n > 3 . (68) 

However it is easy to realize that this cannot be the case, because the theory would not make sense 
as an effective field theory. Following the same logic we used for Ghost Condensation, we find that 
the scaling dimension of the operator tt would be 

_l_l__3_ 

vr ^ s 2"^2n7r . (69) 

This implies that the operator 7r(V7r)^, which is linked by symmetry to the time kinetic term tt^, 
has dimension (7 — 3n)/(2n). For n > 3 this operator is strong at low energy, so that the effective 
field theory does not make sense. 



5 Conclusions 

Given the ongoing experimental effort to test inflation and the proliferation of different models, it is 
quite important to characterize the most general theory of inflation. In this paper we took a novel 
point of view: instead of writing down a general Lagrangian and study the fluctuations around an 
inflating solution, we directly describe the effective theory of fluctuations around a quasi de Sitter 
background, where spatial diffeomorphisms are explicit and the time ones are non-linearly realized. 
We showed that the most generic action can be written at leading order in derivatives in the form 

S = fd-^x ^[Im^iB + A4,Hg«> - A4,{3H^ + H) + ^M2{tf(g'« + l)^ + ^Ah[tf(g'« + 1)3 + 

-^(,» + l)JA-.„ - MlsK^/ - 'MlsK^^SK', + ...] . (70) 

Cosmological correlation functions test this effective field theory at a scale of order the Hubble 
parameter H. In this approach the role of symmetries is made much more transparent. One can 
see explicitly which features are implied by the inflating background solution and in particular the 
quite different behavior in the cases H < 0, H = and H > as the coefficient of the operator g^^ 
is fixed by H. From this point of view, our approach makes clearer the relationship among inflation, 
theories of modification of gravity and theories which violate the Null Energy Condition (equivalent 
to > in the cosmological context) like the bouncing models ^ 1291 [30] . Another example of the 
role of symmetries is given by the link between a reduced speed of sound and an enhanced level of 
non-Gaussianity as both come from the same operator M2{t)^{g^^ + 1)^ and are thus related by the 
non- linear realization of time diffeomorphisms. 
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All the possible deviations from a vanilla slow-roll scenario are systematically encoded in the size 
of higher order operators, similarly to what happens in the study of the Standard Model of particle 
physics. Moreover all single field models are unified in a common framework and this allows us to 
draw general conclusions which are independent of the specific realization, as done in [12\ [3T] for 
example. 

It is easy to think about possible extensions of our formalism. Along the same lines it would 
be interesting to study the most general theory of (single field) quintessence and to work out its 
phenomenological consequences. Differently from inflation, which probes the effective theory at a 
scale of order H, we would be interested in this case to the subhorizon dynamics of perturbations. It 
would also be interesting to use our approach for the study of fluctuations in fluids like in radiation 
or matter dominance [8j. Finally it should be straightforward to introduce additional fields into the 
game and study multi-field infiationary models. 
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Appendix 

A The most general Lagrangian in unitary gauge 

Let us study what are the rules for writing down the most general Lagrangian in unitary gauge. In a 
theory which is only invariant under spatial diffeomorphisms there is a preferred slicing of spacetime 
given by a function t{x) (with time-like gradient), which non- linearly realizes time diffeomorphisms. 
For example if the breaking is given by a time evolving scalar, surfaces of constant t are also of 
constant value of the scalar. Unitary gauge is the one in which the time coordinate t is chosen to 
coincide with t, so that the additional degree of freedom t does not explicitly appear in the action. 
One can therefore build various terms: 

1. Terms which are invariant under all diffeomorphisms: these are just polynomials of the Rie- 
mann tensor R,^,^ and of its covariant derivatives, contracted to give a scalafi 

2. A generic function of t becomes f{t) in unitary gauge. We are therefore free to use generic 
functions of time in front of any terms in the Lagrangian. 

3. The gradient dfj_i becomes 6^ in unitary gauge. Thus in every tensor we can always leave 
free an upper index. For example we can use g^^ (and functions of it) in the unitary gauge 
Lagrangian, or the component of the Ricci tensor R^^. 

^•^The metric and the completely antisymmetric tensor (— g)^i/2gM''po" q^j^ \,q used to contract indices. 
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4. It is useful to define a unit vector perpendicular to surfaces of constant i 



(71) 



This allows to define the induced spatial metric on surfaces of constant t: h^u = g^i, + n^n^. 
Every tensor can be projected on the surfaces using h^^^y. In particular we can use in our action 
the Riemann tensor of the induced 3d metric ^^^Ra/s-yS and covariant derivatives with respect 
to the 3d metric. 

5. Additional possibilities will come from the covariant derivatives of Notice that we can 
equivalently look at covariant derivatives of n^: the derivative acting on the normalization 
factor just gives terms like d^g^^ which are covariant on their own and can be used in the 
unitary gauge Lagrangian. The covariant derivative of projected on the surfaces of constant 
t gives the extrinsic curvature of these surfaces 

K^, = hlV^n, . (72) 

The index v is already projected on the surface because n'^ V^riu = |V^(n^n^) = 0. The 
covariant derivative of perpendicular to the surface can be rewritten as 

n'^V.n, = -l(V°)"'/i(;5^(-5°°) (73) 

so that it does not give rise to new terms. Therefore all covariant derivatives of can be 
written using the extrinsic curvature K^^^ (and its covariant derivatives) and derivatives of 

6. Notice that using at the same time the Riemann tensor of the induced 3d metric and the 
extrinsic curvature is redundant as ^^^Ra/s-yS can be rewritten using the Gauss-Codazzi relation 
as [32] 

^^^Rap-yS = h^h''phPh''sR^,,pa - K^-yKps + Kp^Ko^s ■ (74) 

Thus one can forget about the 3d Riemann tensor altogether. We can also avoid using the 
induced metric h^/s explicitly: written in terms of the 4d metric and one gets only terms 
already discussed abovj^. Finally also the use of covariant derivatives with respect to the 
induced 3d metric can be avoided: the 3d covariant derivative of a projected tensor can be 
obtained as the projection of the 4d covariant derivative [32]. 



We conclude that the most generic action in unitary gauge is given by 

S = J d^xV^ F{Rp,p^,g''°,K^,,V^, t) , (75) 
where all the free indices inside the function F must be upper O's. 



^''Notice that the determinant of the induced metric is related to the one of the full 4d metric hy h — g^^ ■ g 
and that the completely antisymmetric 3d tensor can be rewritten in terms of the 4d one as /i^^/^e*^*^ = 
(-.9)"^^^(-.9°")"^^^e°'^''- 
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B Expanding around a given FRW solution 



In this Section we want to prove that the most generic theory with broken time diffeomorphisms 
around a given FRW background (with k = — 1,0, +1 depending of the spatial curvature) can be 
written as 



(76) 



where the dots stand for terms which are invariant under spatial diffeomorphisms and of quadratic 
(or higher) order in the fluctuations around the given FRW backgrouncJ^. 

As such this statement is trivial. We know that the displayed terms give rise to the wanted 
FRW evolution so that, if we do not want to move away from it, the additional operators must start 
quadratic around this solution. What we want to say is that each one of the additional invariant 
terms is quadratic (or of higher order) without cancellation of linear contributions among various 
operators. These terms will be written as polynomials (quadratic and higher) of linear operators 
like g^^ + 1, 5Kfj_^ = Kfj_^ — K^^} ■, ^R^ivpa = Rpupa — R^pupa and so on. Notice that these terms start 
linear in the perturbations as we have explicitly removed their value evaluated on the given FRW 
solution. Given the symmetries of a FRW metric, every tensor evaluated on the background {K^} , 
Rp^upa, i^aRpi^pa)^^^ •••) cau be Written just in terms of g^u, f^p and functions of time. For example 

4°) = a'Hh^, (77) 
= 2{H + k)h^^pK^, + {H + H^)a^h^jyp + V^^rn. (78) 

where A; is a constant which depends on the curvature of the spatial slices and the permutations 
are acting only on the last term. As such all the operators evaluated on the FRW background are 
themselves covariant operators, so that operators like 5Kfj_y and 5Rf^upa are well defined covariant 
operators which vanish on the given FRW background and start linear in the perturbations. We 
stress that this possibility of rewriting the tensors evaluated on the background holds only because of 
the high degree of symmetry of the FRW background and it would not be true if one were interested 
in expanding around a less symmetric solution, e.g. a non-homogeneous background. 

In equation (|76p only the displayed operators contain linear terms in the fluctuations, so that 
the coefficients of \/—g g^^ and \/—g are uniquely determined by the background FRW solution. 

Let us now see how the Lagrangian can always be cast in the form (I76p . If we take an opera- 
tor composed by the contraction of two tensors T and G (the generalization with more tensors is 
straightforward) we can write 

TG = 6T5G + T^^^G + TG^^^ - T^^^G^^^ . (79) 

Let us discuss each term of the sum. The first one starts explicitly quadratic in the perturbation 
as we want. As we said, given the symmetries of the FRW background, the unperturbed tensors 
r(0) and can be written as functions of g^^p., and t. Therefore the last term T^'^'^G^^'^ is just 
a polynomial of g^^ with time dependent coefficients; it contains the terms ^J—g g^^ and \/—g plus 
operators which start explicitly quadratic in the perturbations. We are left with tensors of the form 



^^We can always make the coefficient in front of R time independent through an appropriate field redefinition 
9pL' ~^ 9fiL' ■ f{t)- This corresponds, in the usual formalism, to going to Einstein frame. 
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T^^^G. We want to prove that also these terms can be written as the linear operators in eq. (j76p 
plus operators that start quadratic in the fluctuations. By construction G will be linear either in 
K/^ii, or Rj_ippa with covariant derivatives acting on them. Covariant derivatives can be dealt with 
by successive integration by parts, letting them act on T^'^^ and the time dependent coefficient of 
the operator. In doing so we can generate extrinsic curvature terms. In this case we can reiterate 
eq. (|79p until no covariant derivatives are lefj^. We are thus left with the only possible scalar linear 
terms with no covariant derivatives: ^ and R^^ . Both of them can be rewritten in a more useful 
form. We can integrate by parts the extrinsic curvature term 

j d^x^ f{t)K^^ = J d^x^ f V^n^ = - j d^x^ nf^d^f = J d^x^^/^f . (80) 
While we can deal with R^^ using the following relationship 



(-5°°)-^i?°° = R^^un^n'' = K^- K^.K^^'^ - V^(n'^V,n'^) + V,(n'^V^n'') . (81) 
The last two terms can again be integrated by parts: 

j d^x^fm^in'^V.nn = - j d^x^ d^f n^K% , (82) 



j d^x^f{t)V,{n^'V^n'') = - j d^x^ d,f n^'V^n'^ = 



(83) 



where in the last passage we have used that d^f oc n^. This shows that Kj^ and R^^ can be written 
in terms of the linear operators of eq. (|76|) plus invariant terms that starts quadratically in the 
fluctuations. 

In conclusion, we have shown that the most general Lagrangian of a theory with broken time 
difFeomorphisms around a given FRW background can be written in the form: 



5 



d^xy/—g 



1 



k 



-Mi^R + M^AH - ^) ■ g^^ - MMH' + H + 2^] + 



k 



(84) 



F(2)(500^1,<5i^^,,<5i?^,^,;V^;t) 
where F^^^ starts quadratic in the arguments g'^^ + 1, SK^i, and SR^j^ppa- 
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